Viper Venom Phospholipase A2 Database: The Structural and Functional Anatomy of a Primary Toxin in Envenomation.

Viper venom phospholipase A2 enzymes (vvPLA2s) and phospholipase A2-like (PLA2-like) proteins are two of the principal toxins in viper venom that are responsible for the severe myotoxic and neurotoxic effects caused by snakebite envenoming, among other pathologies. As snakebite envenoming is the deadliest neglected tropical disease, a complete understanding of these proteins' properties and their mechanisms of action is urgently needed. Therefore, we created a database comprising information on the holo-form, cofactor-bound 3D structure of 217 vvPLA2 and PLA2-like proteins in their physiologic environment, as well as 79 membrane-bound viper species from 24 genera, which we have made available to the scientific community to accelerate the development of new anti-snakebite drugs. In addition, the analysis of the sequenced, 3D structure of the database proteins reveals essential aspects of the anatomy of the proteins, their toxicity mechanisms, and the conserved binding site areas that may anchor universal interspecific inhibitors. Moreover, it pinpoints hypotheses for the molecular origin of the myotoxicity of the PLA2-like proteins. Altogether, this study provides an understanding of the diversity of these toxins and how they are conserved, and it indicates how to develop broad, interspecies, efficient small-molecule inhibitors to target the toxin's many mechanisms of action.

Deciphering the Catalytic Mechanism of Virginiamycin B Lyase with Multiscale Methods and Molecular Dynamics Simulations.

Due to the emergence of antibiotic resistance, the need to explore novel antibiotics and/or novel strategies to counter antibiotic resistance is of utmost importance. In this work, we explored the molecular and mechanistic details of the degradation of a streptogramin B antibiotic by virginiamycin B (Vgb) lyase of Staphylococcus aureus using classical molecular dynamics simulations and multiscale quantum mechanics/molecular mechanics methods. Our results were in line with available experimental kinetic information. Although we were able to identify a stepwise mechanism, in the wild-type enzyme, the intermediate is short-lived, showing a small barrier to decay to the product state. The impact of point mutations on the reaction was also assessed, showing not only the importance of active site residues to the reaction catalyzed by Vgb lyase but also of near positive and negative residues surrounding the active site. Using molecular dynamics simulations, we also predicted the most likely protonation state of the 3-hydroxypicolinic moiety of the antibiotic and the impact of mutants on antibiotic binding. All this information will expand our understanding of linearization reactions of cyclic antibiotics, which are crucial for the development of novel strategies that aim to tackle antibiotic resistance.

Unraveling the Reaction Mechanism of Russell's Viper Venom Factor X Activator: A Paradigm for the Reactivity of Zinc Metalloproteinases?

Snake venom metalloproteinases (SVMPs) are important drug targets against snakebite envenoming, the neglected tropical disease with the highest mortality worldwide. Here, we focus on Russell's viper (Daboia russelii), one of the "big four" snakes of the Indian subcontinent that, together, are responsible for ca. 50,000 fatalities annually. The "Russell's viper venom factor X activator" (RVV-X), a highly toxic metalloproteinase, activates the blood coagulation factor X (FX), leading to the prey's abnormal blood clotting and death. Given its tremendous public health impact, the WHO recognized an urgent need to develop efficient, heat-stable, and affordable-for-all small-molecule inhibitors, for which a deep understanding of the mechanisms of action of snake's principal toxins is fundamental. In this study, we determine the catalytic mechanism of RVV-X by using a density functional theory/molecular mechanics (DFT:MM) methodology to calculate its free energy profile. The results showed that the catalytic process takes place via two steps. The first step involves a nucleophilic attack by an in situ generated hydroxide ion on the substrate carbonyl, yielding an activation barrier of 17.7 kcal·mol-1, while the second step corresponds to protonation of the peptide nitrogen and peptide bond cleavage with an energy barrier of 23.1 kcal·mol-1. Our study shows a unique role played by Zn2+ in catalysis by lowering the pKa of the Zn2+-bound water molecule, enough to permit the swift formation of the hydroxide nucleophile through barrierless deprotonation by the formally much less basic Glu140. Without the Zn2+ cofactor, this step would be rate-limiting.

Beyond the TPP+ "gold standard": a new generation mitochondrial delivery vector based on extended PN frameworks.

Mitochondrial targeting represents an attractive strategy for treating metabolic, degenerative and hyperproliferative diseases, since this organelle plays key roles in essential cellular functions. Triphenylphosphonium (TPP+) moieties - the current "gold standard" - have been widely used as mitochondrial targeting vectors for a wide range of molecular cargo. Recently, further optimisation of the TPP+ platform drew considerable interest as a way to enhance mitochondrial therapies. However, although the modification of this system appears promising, the core structure of the TPP+ moiety remains largely unchanged. Thus, this study explored the use of aminophosphonium (PN+) and phosphazenylphosphonium (PPN+) main group frameworks as novel mitochondrial delivery vectors. The PPN+ moiety was found to be a highly promising platform for this purpose, owing to its unique electronic properties and high lipophilicity. This has been demonstrated by the high mitochondrial accumulation of a PPN+-conjugated fluorophore relative to its TPP+-conjugated counterpart, and has been further supported by density functional theory and molecular dynamics calculations, highlighting the PPN+ moiety's unusual electronic properties. These results demonstrate the potential of novel phosphorus-nitrogen based frameworks as highly effective mitochondrial delivery vectors over traditional TPP+ vectors.

Seed Size, Not Dispersal Syndrome, Determines Potential for Spread of Ricefield Weeds by Gulls.

Recent field data suggest that migratory gulls disperse many rice field weeds by gut passage (endozoochory), most of which are dry fruited and widely assumed to have no long-distance dispersal mechanisms, except via human activity. We investigated this mechanism with a feeding experiment, in which seeds of five common rice field weeds (in order of increasing seed size: Juncus bufonius, Cyperus difformis, Polypogon monspeliensis, Amaranthus retroflexus, and the fleshy-fruited Solanum nigrum) were fed to seven individuals of lesser black-backed gulls Larus fuscus held in captivity. We quantified seed survival after collecting faeces at intervals for 33 h after ingestion, then extracting intact seeds and running germination tests, which were also conducted for control seeds. All five species showed high seed survival after gut passage, of >70%. Gut retention times averaged 2-4 h, but maxima exceeded 23 h for all species. Germinability after gut passage was 16-54%, and gut passage accelerated germination in J. bufonius and S. nigrum, but slowed it down in the other species. All species had lower germinability after gut passage compared to control seeds (likely due to stratification prior to the experiment), but the loss of germinability was higher in smaller seeds. There was no evidence that the different dispersal syndromes assigned to the five species (endozoochory, epizoochory or barochory) had any influence on our results. In contrast, mean gut retention time was strongly and positively related to seed size, likely because small seeds pass more quickly from the gizzard into the intestines. Non-classical endozoochory of dry-fruited seeds by waterbirds is a major but overlooked mechanism for potential long-distance dispersal, and more research into this process is likely essential for effective weed management.

Development of Nanoscale Graphene Oxide Models for the Adsorption of Biological Molecules.

Graphene oxide (GO), a nanomaterial with promising applications that range from water purification to enzyme immobilization, is actively present in scientific research since its discovery. GO studies with computational methodologies such as molecular dynamics are frequently reported in the literature; however, the models used often rely on approximations, such as randomly placing functional groups and the use of generalized force fields. Therefore, it is important to develop new MD models that provide a more accurate description of GO structures and their interaction with an aqueous solvent and other adsorbate molecules. In this paper, we derived new force field non-bonded parameters from linear-scaling density functional theory calculations of nanoscale GO sheets with more than 10,000 atoms through an atoms-in-molecules (AIM) partitioning scheme. The resulting GAFF2-AIM force field, derived from the bonded terms of GAFF2 parameterization, reproduces the solvent structure reported in ab initio MD simulations better than the force field nowadays widely used in the literature. Additionally, we analyzed the effect of the ionic strength of the medium and of the C/O ratio on the distribution of charges surrounding the GO sheets. Finally, we simulated the adsorption of natural amino acid molecules to a GO sheet and estimated their free energy of binding, which compared very favorably to their respective experimental values, validating the force field presented in this work.

Engineering DszC Mutants from Transition State Macrodipole Considerations and Evolutionary Sequence Analysis.

We describe an approach to identify enzyme mutants with increased turnover using the enzyme DszC as a case study. Our approach is based on recalculating the barriers of alanine mutants through single-point energy calculations at the hybrid QM/MM level in the wild-type reactant and transition state geometries. We analyze the difference in the electron density between the reactant and transition state to identify sites/residues where electrostatic interactions stabilize the transition state over the reactants. We also assess the insertion of a unit probe charge to identify positions in which the introduction of charged residues lowers the barrier.

QM/MM Study of the Reaction Mechanism of Thermophilic Glucuronoyl Esterase for Biomass Treatment.

Hydrolysis of lignocellulosic biomass, composed of a lignin-carbohydrate-complex (LCC) matrix, is critical for producing bioethanol from glucose. However, current methods for LCC processing require costly and polluting processes. The fungal Thermothelomyces thermophila glucuronoyl esterase (TtGE) is a promising thermophilic enzyme that hydrolyses LCC ester bonds. This study describes the TtGE catalytic mechanism using QM/MM methods. Two nearly-degenerate rate-determining transition states were found, with barriers of 16 and 17 kcal ⋅ mol-1 , both with a zwitterionic nature that results from a proton interplay from His346 to either the Ser213-hydroxyl or the lignin leaving group and the rehybridisation of the ester moiety of the substrate to an alkoxide. An oxyanion hole, characteristic of esterases, was provided by the conserved Arg214 through its backbone and sidechain. Our work further suggests that a mutation of Glu267 to a non-negative residue will decrease the energetic barrier in ca. -5 kcal ⋅ mol-1 , improving the catalytic rate of TtGE.

Role of Enzyme and Active Site Conformational Dynamics in the Catalysis by α-Amylase Explored with QM/MM Molecular Dynamics.

We assessed enzyme:substrate conformational dynamics and the rate-limiting glycosylation step of a human pancreatic α-amylase:maltopentose complex. Microsecond molecular dynamics simulations suggested that the distance of the catalytic Asp197 nucleophile to the anomeric carbon of the buried glucoside is responsible for most of the enzyme active site fluctuations and that both Asp197 and Asp300 interact the most with the buried glucoside unit. The buried glucoside binds either in a 4C1 chair or 2SO skew conformations, both of which can change to TS-like conformations characteristic of retaining glucosidases. Starting from four distinct enzyme:substrate complexes, umbrella sampling quantum mechanics/molecular mechanics simulations (converged within less than 1 kcal·mol-1 within a total simulation time of 1.6 ns) indicated that the reaction occurrs with a Gibbs barrier of 13.9 kcal·mol -1, in one asynchronous concerted step encompassing an acid-base reaction with Glu233 followed by a loose SN2-like nucleophilic substitution by the Asp197. The transition state is characterized by a 2H3 half-chair conformation of the buried glucoside that quickly changes to the E3 envelope conformation preceding the attack of the anomeric carbon by the Asp197 nucleophile. Thermodynamic analysis of the reaction supported that a water molecule tightly hydrogen bonded to the glycosidic oxygen of the substrate at the reactant state (∼1.6 Å) forms a short hydrogen bond with Glu233 at the transition state (∼1.7 Å) and lowers the Gibbs barrier in over 5 kcal·mol-1. The resulting Asp197-glycosyl was mostly found in the 4C1 conformation, although the more endergonic B3,O conformation was also observed. Altogether, the combination of short distances for the acid-base reaction with the Glu233 and for the nucleophilic attack by the Asp197 nucleophile and the availability of water within hydrogen bonding distance of the glycosidic oxygen provides a reliable criteria to identify reactive conformations of α-amylase complexes.

The chemistry of snake venom and its medicinal potential.

The fascination and fear of snakes dates back to time immemorial, with the first scientific treatise on snakebite envenoming, the Brooklyn Medical Papyrus, dating from ancient Egypt. Owing to their lethality, snakes have often been associated with images of perfidy, treachery and death. However, snakes did not always have such negative connotations. The curative capacity of venom has been known since antiquity, also making the snake a symbol of pharmacy and medicine. Today, there is renewed interest in pursuing snake-venom-based therapies. This Review focuses on the chemistry of snake venom and the potential for venom to be exploited for medicinal purposes in the development of drugs. The mixture of toxins that constitute snake venom is examined, focusing on the molecular structure, chemical reactivity and target recognition of the most bioactive toxins, from which bioactive drugs might be developed. The design and working mechanisms of snake-venom-derived drugs are illustrated, and the strategies by which toxins are transformed into therapeutics are analysed. Finally, the challenges in realizing the immense curative potential of snake venom are discussed, and chemical strategies by which a plethora of new drugs could be derived from snake venom are proposed.

Different Enzyme Conformations Induce Different Mechanistic Traits in HIV-1 Protease.

The influence of the dynamical flexibility of enzymes on reaction mechanisms is a cornerstone in biological sciences. In this study, we aim to 1) study the convergence of the activation free energy by using the first step of the reaction catalysed by HIV-1 protease as a case study, and 2) provide further evidence for a mechanistic divergence in this enzyme, as two different reaction pathways were seen to contribute to this step. We used quantum mechanics/molecular mechanics molecular dynamics simulations, on four different initial conformations that led to different barriers in a previous study. Despite the sampling, the four activation free energies still spanned a range of 5.0 kcal ⋅ mol-1 . Furthermore, the new simulations did confirm the occurrence of an unusual mechanistic divergence, with two different mechanistic pathways displaying equivalent barriers. An active-site water molecule is proposed to influence the mechanistic pathway.

Towards the Accurate Thermodynamic Characterization of Enzyme Reaction Mechanisms.

We employed QM/MM molecular dynamics (MD) simulations to characterize the rate-limiting step of the glycosylation reaction of pancreatic α-amylase with combined DFT/molecular dynamics methods (PBE/def2-SVP : AMBER). Upon careful choice of four starting active site conformations based on thorough reactivity criteria, Gibbs energy profiles were calculated with umbrella sampling simulations within a statistical convergence of 1-2 kcal ⋅ mol-1 . Nevertheless, Gibbs activation barriers and reaction energies still varied from 11.0 to 16.8 kcal ⋅ mol-1 and -6.3 to +3.8 kcal ⋅ mol-1 depending on the starting conformations, showing that despite significant state-of-the-art QM/MM MD sampling (0.5 ns/profile) the result still depends on the starting structure. The results supported the one step dissociative mechanism of Asp197 glycosylation preceded by an acid-base reaction by the Glu233, which are qualitatively similar to those from multi-PES QM/MM studies, and thus support the use of the latter to determine enzyme reaction mechanisms.

Necessity is the Mother of Invention: A Remote Molecular Bioinformatics Practical Course in the COVID-19 Era.

The COVID-19 pandemic has brought many challenges to human beings, related to not only health and way of life but also teaching because of the interruption of the standard training at universities imposed by lockdowns. Concerning the latter, the academic community had to reinvent itself, in many ways, to carry on with prepandemic education. This article focuses on the use of modern technology and software to create a virtual, highly interactive classroom where a remote but still hands-on course on molecular bioinformatics can be taught, motivating the university students and helping them learn the course contents without significant compromises imposed by successive lockdowns. We implemented such a virtual hands-on molecular bioinformatics course in the second semester of the 2020/2021 academic year. Furthermore, we compared the learning outcomes with those for the earlier editions of the same course in the pre-COVID-19 era, in which the more traditional teaching method was used where all teaching was delivered with physically present lecturers. The virtual classroom proposed here allowed the students to develop skills close to, although slightly below, those obtained with physically present learning.

Root-knot nematodes (Meloidogyne spp.) a threat to agriculture in Mexico: biology, current control strategies, and perspectives.

Root-knot nematodes (RKN) are sedentary parasites of the roots of plants and are considered some of the most damaging pests in agriculture. Since RKN target the root vascular system, they provoke host nutrient deprivation and defective water transport, causing above-ground symptoms of growth stunting, wilting, chlorosis, and reduced crop yields. In Mexico RKN infestations are primarily dealt with by treating with synthetic chemically based nematicides that are preferred by farmers over available bioproducts. However, due to environmental and human health concerns chemical control is increasingly restricted. Biological control of RKNs can help reduce the use of chemical nematicides as it is achieved with antagonistic organisms, mainly bacteria, fungi, other nematodes, or consortia of diverse microorganisms, which control nematodes directly by predation and parasitism at different stages: eggs, juveniles, or adults; or indirectly by the action of toxic diffusible inhibitory metabolites. The need to increase agricultural production and reduce negative environmental impact creates an opportunity for optimizing biological control agents to suppress nematode populations, but this endeavour remains challenging as researchers around the world try to understand diverse control mechanisms, nematode and microbe life cycles, ecology, metabolite production, predatory behaviours, molecular and biochemical interactions, in order to generate attractive products with the approval of local regulatory bodies. Here, we provide a brief review of the biology of the genus Meloidogyne, biological control strategies, and a comparison between chemical and bioproducts in the Mexican market, and guidelines emitted by national agencies to ensure safety and effectiveness of new developments.

The chemistry of snake venom and its medicinal potential.

The fascination and fear of snakes dates back to time immemorial, with the first scientific treatise on snakebite envenoming, the Brooklyn Medical Papyrus, dating from ancient Egypt. Owing to their lethality, snakes have often been associated with images of perfidy, treachery and death. However, snakes did not always have such negative connotations. The curative capacity of venom has been known since antiquity, also making the snake a symbol of pharmacy and medicine. Today, there is renewed interest in pursuing snake-venom-based therapies. This Review focuses on the chemistry of snake venom and the potential for venom to be exploited for medicinal purposes in the development of drugs. The mixture of toxins that constitute snake venom is examined, focusing on the molecular structure, chemical reactivity and target recognition of the most bioactive toxins, from which bioactive drugs might be developed. The design and working mechanisms of snake-venom-derived drugs are illustrated, and the strategies by which toxins are transformed into therapeutics are analysed. Finally, the challenges in realizing the immense curative potential of snake venom are discussed, and chemical strategies by which a plethora of new drugs could be derived from snake venom are proposed.

Modern Strategies for the Diversification of the Supply of Natural Compounds: The Case of Alkaloid Painkillers.

Plant-derived natural compounds have been used for treating diseases since prehistorical times. The supply of many plant-derived natural compounds for medicinal purposes, such as thebaine, morphine, and codeine, is primarily dependent on opium poppy crop harvesting. This dependency adds an extra risk factor to ensuring the supply chain because crops are highly susceptible to environmental conditions. Emerging technologies, such as biocatalysis, might help to solve this problem by diversifying the sources of supply of these compounds. Here we review the first committed step in the production of alkaloid painkillers, the production of S-norcoclaurine, and the enzymes involved. The improvement of these enzymes can be carried out experimentally by directed evolution and rational design strategies, supported by computational methods, to create variants that produce the S-norcoclaurine precursor for alkaloid painkillers in heterologous organisms, meeting the pharmaceutical industry standards and needs without depending on opium poppy crops.

Exploring the permeation of fluoroquinolone metalloantibiotics across outer membrane porins by combining molecular dynamics simulations and a porin-mimetic in vitro model.

The misuse and overuse of fluoroquinolones in recent years have triggered alarming levels of resistance to these antibiotics. Porin channels are crucial for the permeation of fluoroquinolones across the outer membrane of Gram-negative bacteria and modifications in porin expression are an important mechanism of bacterial resistance. One possible strategy to overcome this problem is the development of ternary copper complexes with fluoroquinolones. Compared to fluoroquinolones, these metalloantibiotics present a larger partition to the lipid bilayer and a more favorable permeation, by passive diffusion, across bacteriomimetic phospholipid-based model membranes. To rule out the porin-dependent pathway for the metalloantibiotics, we explored the permeation through OmpF (one of the most abundant porins present in the outer membrane of Gram-negative bacteria) using a multi-component approach. X-ray studies of OmpF porin crystals soaked with a ciprofloxacin ternary copper complex did not show a well-defined binding site for the compound. Molecular dynamics simulations showed that the translocation of the metalloantibiotic through this porin is less favorable than that of free fluoroquinolone, as it presented a much larger free energy barrier to cross the narrow constriction region of the pore. Lastly, permeability studies of different fluoroquinolones and their respective copper complexes using a porin-mimetic in vitro model corroborated the lower rate of permeation for the metalloantibiotics relative to the free antibiotics. Our results support a porin-independent mechanism for the influx of the metalloantibiotics into the bacterial cell. This finding brings additional support to the potential application of these metalloantibiotics in the fight against resistant infections and as an alternative to fluoroquinolones.

Structure based virtual screening of natural product molecules as glycosidase inhibitors.

Objective of the present investigation comprised of the application of in silico methods to discover novel natural product (NP) based potential inhibitors for carbohydrate mediated diseases. Structure based drug design studies (molecular docking and structure based pharmacophore analysis)  were carried out on a series of natural product compounds to identify significant bioactive molecules to inhibit α-mannosidase (I and II) and ß-galactosidase enzymes. Furthermore, protein ligand interaction fingerprint analysis, molecular dynamics simulations and molecular access system (MACCS) fingerprint analysis were performed to understand the binding behaviors of the studied molecules. The results derived from these analyses showed that the identified compounds exhibit significant binding interactions with the active site residues. The compounds, NP-51, NP-81 and NP-165 have shown significant docking score against the studied enzymes (α-mannosidases-I, α-mannosidases-II and ß-galactosidases). The fingerprint studies showed that the presence of rings (aromatic or aliphatic) with sulfur atoms, nitrogen atoms, methyl groups, etc. have favorable effects on the α-mannosidase II inhibitory activity. However, the presence of halogen atoms substituted in the molecules have reduced inhibitory ability against α-mannosidase II. The compound, NP-165 has significant activity against both enzymes (α-mannosidases and ß-galactosidases). These studies accomplished that the compounds identified through in silico methodologies can be used to develop semisynthetic derivatives of the glycosidase inhibitors and can be screened for the treatment of different carbohydrate mediated diseases. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s40203-021-00115-9.

The Catalytic Mechanism of the Retaining Glycosyltransferase Mannosylglycerate Synthase.

To protect their intracellular proteins, extremophile microorganisms synthesize molecules called compatible solutes. These molecules are the result of the attachment of a small negatively charged molecule to a sugar molecule. It has been found that these molecules, not only protect the microorganism against osmotic stress but also against other extreme conditions. They can also confer protection against extreme conditions to isolated enzymes from different organisms making them an exciting prospect for potential biotechnological applications. One of the most widespread compatible solute in hyperthermophile organisms is the molecule 2-O-α-D-mannosyl-D-glycerate (MG). In addition to confer protection to proteins against extreme conditions, MG was found to prevent Alzheimer's ß-amyloid aggregation and reduce α-synuclein fibril formation in Parkinson's disease. In this work we studied, using computational methods, the catalytic mechanism of the synthesis of MG by the enzyme mannosylglycerate synthase (MGS) from the thermophilic bacteria Rhodothermus marinus.

Animal Fatty Acid Synthase: A Chemical Nanofactory.

Fatty acids are crucial molecules for most living beings, very well spread and conserved across species. These molecules play a role in energy storage, cell membrane architecture, and cell signaling, the latter through their derivative metabolites. De novo synthesis of fatty acids is a complex chemical process that can be achieved either by a metabolic pathway built by a sequence of individual enzymes, such as in most bacteria, or by a single, large multi-enzyme, which incorporates all the chemical capabilities of the metabolic pathway, such as in animals and fungi, and in some bacteria. Here we focus on the multi-enzymes, specifically in the animal fatty acid synthase (FAS). We start by providing a historical overview of this vast field of research. We follow by describing the extraordinary architecture of animal FAS, a homodimeric multi-enzyme with seven different active sites per dimer, including a carrier protein that carries the intermediates from one active site to the next. We then delve into this multi-enzyme's detailed chemistry and critically discuss the current knowledge on the chemical mechanism of each of the steps necessary to synthesize a single fatty acid molecule with atomic detail. In line with this, we discuss the potential and achieved FAS applications in biotechnology, as biosynthetic machines, and compare them with their homologous polyketide synthases, which are also finding wide applications in the same field. Finally, we discuss some open questions on the architecture of FAS, such as their peculiar substrate-shuttling arm, and describe possible reasons for the emergence of large megasynthases during evolution, questions that have fascinated biochemists from long ago but are still far from answered and understood.